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5— ( ! ABSTRACT 



< 



Profile comparison of the Stokes parameters V and I is a powerful tool for maser data 
analysis, providing the first direct methods for unambiguous determination of (1) the 
maser saturation stage, (2) the amplification optical depth and intrinsic Doppler width 
of unsaturated masers, and (3) the comparative magnitudes of Zeeman splitting and 
Doppler linewidth. Circular polarization recently detected in OH 1720 MHz emission 
from the Galactic center appears to provide the first direct evidence for maser saturation. 
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INTRODUCTION 
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■ Few properties of astronomical masers are determined directly by observations, most are in- 

ferred indirectly. Foremost among the latter is the maser saturation stage. Saturation has a 
significant impact on maser growth, so determining whether a maser is saturated (J > J s , where J 
is the angle-averaged intensity and J s is the saturation intensity) or not is usually a precondition 
for analysis of the observations. Unfortunately, this crucial issue is not convincingly settled. Strong 
masers are generally believed to be saturated, but in most cases the evidence is less than compelling 
as it relies primarily on plausibility arguments rather than quantitative tests (see e.g. Elitzur 1992| , 



§8.6). This unsatisfactory situation reflects a fundamental difficulty — neither J nor J s is directly 
measurable. The saturation parameter J s is a theoretical quantity, determined only within the 
context of a given pumping scheme. And because maser radiation is highly beamed, J = I£l/Air 
so this quantity, too, cannot be directly measured; the intensity / is measurable when the maser is 
resolved, but the beaming angle f2 is unobservable. Similarly, the amplification optical depth has 
never been directly determined for any maser that amplifies its own radiation. 



Recent VLA observations of OH 1720 MHz masers near the Galactic center by Yusef-Zadeh 



et al. 1996| open up new possibilities for direct determination of some maser properties. Significant 



circular polarization (upward of 20%) is detected in various spectral features, and the right- and 
left-hand components coincide on the sky, as expected from the Zeeman effect. Furthermore, the 
spectral shape of the Stokes parameter V follows an antisymmetric S-curve with sharp reversal 
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at line center, the typical profile for Zeeman shift Aub much smaller than the Doppler linewidth 



Avd- Similar results were previously reported for H2O masers in star-forming regions by Fiebig & 



Gusten 1989 and for OH 1612 MHz masers in OH/IR stars by Zell & Fix 1991, but the polarization 
was lower and the quality of the data not nearly as high. 

The general maser polarization solution was recently derived for arbitrary values of 

Au B , , 



(Elitzur 1996, hereafter E96) and the solution properties at < 1 closely match the observed 



circular polarization. Here I show that a comparative analysis of the spectral profiles of I and V, two 
measurable independent maser intensities, offers direct determination of various maser properties, 
in particular the saturation stage. The analysis is readily performed with the aid of the ratio profile 

K(v) - - V ^ (2) 

where the prime denotes derivative with respect to frequency and v = V/I is the fractional circular 
polarization. When < 1, spectral analysis of 1Z offers intrinsic sensitivity of order xb and 
has long been an important tool in studies of the Zeeman effect of thermal radiation (see e.g. 
Troland fc Heiles 1982j ). In that case 1Z is constant across the spectral line and its magnitude 



determines the magnetic field along the line of sight. This constancy of 1Z follows from some simple, 
general symmetry arguments as shown by prutcher et al. 1993; (see also §3 below). However, maser 
exponential amplification during unsaturated growth destroys both the underlying symmetry and 
the constancy of 1Z, the saturation process restores both. The key to the different behavior, and 
7£-profiles, in the two regimes is the narrowing of the maser line during unsaturated amplification 
and its rebroadening during saturation. 

The important differences between thermal and maser polarization are discussed in detail 



below. For completeness, some basic elements of the polarization theory developed in [E96| are 
reproduced in §2. The 1Z profile is discussed in §3 for thermal radiation and in §4 for maser 
radiation when xb *C 1. In §5, circular polarization for fully resolved Zeeman patterns, xb > 1, is 
discussed. The implications for observations are discussed in detail in §6. 



GENERALITIES 



A full description of electromagnetic radiation involves the 4-vector of its Stokes parameters 
S = (I, Q, U, V). The general transfer equation for S is 



dS_ 
~dl 



e + K-S, 



(3) 



where e is the 4-vector analog of the familiar volume emission coefficient and K is the matrix 
analog of the absorption coefficient. In the case of line radiation, K was derived in Litvak 1975| and 
is reproduced here using the notation of E96. For any spin transition, the radiative interactions 
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with photons of polarization Am (= +1, 0, —1) are characterized by an absorption coefficient n Am . 
Introduce 

k 1 = !(«+ + «-), ( 4 ) 
the mean absorption coefficient for Am = ±1 transitions, and the three linear combinations 

Km = i[ K 1 (l + cos 2 ^) + K O sin 2 0] 
«, = \{k 1 - K °)sin 2 (9 

K c = t;(k + — K~) COS 6, (5) 

where 9 is the angle of the radiation propagation direction from the quantization axis, taken here 
along the magnetic field. Then the transfer matrix for line radiation is 
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(6) 



The polarization of the 4- vector S is characterized by II = (q, u, v), the 3- vector of its normal- 
ized Stokes parameters q = Q/I,u = U/I,v = V/I. It is easy to verify that the transfer matrix K 
preserves full polarization, |II| = 1. From Maxwell's equations, radiation is always emitted fully 
polarized and for any intensity / there is an ensemble of configurations of S (polarization modes) 
that differ from each other only by the directions of their corresponding unit vectors II, i.e., the 
sense of polarization. For each mode, II remains a unit vector and radiative transfer only rotates 
it, transferring polarization between the linear and circular polarizations of individual modes. The 
polarization vector of the entire radiation field at intensity I is obtained from the ensemble average 
of all polarization modes. This average differs from zero, and the observed radiation is polarized, 
only when certain modes are favored. 

The emission term for the mode with initial polarization IIo can be written as 

e = K m S , where S = S (l, U ) (7) 

and Sq is the standard source function. The ensemble average of these source terms has the common 
intensity Sq, and when the individual polarizations point at random directions in q,u,v space it 
also has zero polarization. This is usually the case for radiation generated in spontaneous decays. 
The radiative transfer equations admit the formal solution for the intensity of individual modes 

I = S e Tm+T «, (8) 

where 

The optical depth t w describes the rotation of the mode polarization while the growth of its intensity 
is dominated by r m , which is always larger than r n . Modes that have r n = do not rotate, providing 
a stationary polarization configuration. 
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3. THERMAL CIRCULAR POLARIZATION 

In the case of thermal emission, the emerging radiation samples no more than ~ one optical 
depth into the source and its properties primarily reflect those of the source terms. Therefore, 
only the polarization of the emission terms need be considered; the subsequent effect of radiative 
transfer is generally only a perturbation. 

In the presence of a magnetic field, thermal emission is circularly polarized and the polarization 



can be derived from rather general considerations (e.g., Crutcher et al. 1993). The energies of 
magnetic sub- levels are shifted by the Zeeman effect, splitting each line into components with 
Am = —1, 0, +1 (where Am refers to the m-change in absorption for compatibility with the previous 
section). The standard Zeeman emission pattern is comprised of three components centered on the 
frequencies vq — AmAi/g, each one with appropriate polarization properties. From symmetry, for 
m-independent line excitations the intensities of the three components are related via 

I ± (x)=I (x±x B )~I°(x)±x B ^^- (10) 

ax 

where x = {y — v§)j Avd is the dimensionless frequency shift from line center. The last approxi- 
mation is valid for all x when ib < 1. Because of the incoherence of thermal emission, the Stokes 
parameters are simple superpositions of the component intensities 



I(x) = +r){l + cos 2 #) + 2/°sin 2 6>] ~ I°(x) 

dl°(x) 

V(x) = U I+ ~ r )cos6 ^ x B — 7^COS0 (11) 
2 dx 

where the final expressions are the terms to leading order in xg. Therefore, 

V{v) = Av B caB0I'(v), i.e. H(v) = Av B cos 6 (12) 

(note that x B d/dx = Au B d/ *dv). The constancy of the ratio profile 1Z reflects the symmetry of 
equation [H] — the three I Am are described by a single spectral profile at appropriately shifted 
arguments. This symmetry is unique to thermal radiation. The general symmetry of the problem 
is invariance under mirror reflections perpendicular to the magnetic axis and its consequence is 
instead 

I + {u Q + 8v)=r{u Q -8v) (13) 

for any bv. Equation |l^ follows only if the spectral shape of each I Am is additionally symmetric 
about its centroid at vq — AmAi/g. This is the case for thermal radiation because the emission is 
invariant under reversal of the particle motions. 



4. MASER CIRCULAR POLARIZATION 

In contrast with thermal emission, the Stokes parameters of maser radiation are dominated by 
the interaction terms in eq. [3| rather than the source terms when J ^ \J SqJ s ( E96| ) . This condition 
is obeyed by virtually all bright maser sources. The polarization properties of radiation generated 
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in spontaneous decays and in stimulated emissions are entirely different from each other. Consider 
for example a photon produced in a Am = decay and thus linearly polarized. This photon 
contributes only to the linear polarization of the source terms and will be reflected as such in the 
polarization properties of observed thermal radiation. However, this photon can induce transitions 
with any value of Am, thus any polarization, because it can also be described as a coherent mixture 
of two circularly polarized photons. Consider, for instance, the stimulated emissions induced by 
this photon in the case of a spin 1 — » transition. When the interacting particle is in the m = 
state, the induced photon is linearly polarized because it, too, is generated in a Am = transition. 
But when the particle is in one of the |m| = 1 states, the induced photon is circularly polarized 
because it is generated in a Am = ±1 transition (the interaction amplitude is reduced, though). 
A photon generated in stimulated emission has the same wave vector as the parent photon but 
not necessarily the same polarization, because the latter depends also on the magnetic quantum 
number of the interacting particle. Because of the finite line widths, different Am transitions 
overlap when 15 C 1 and particles in the same magnetic subdevel can interact with photons of 
different polarizations. Amplification by stimulated emission mixes the original polarizations of all 
the photons created in spontaneous decays. 

The polarization mixing effect of the amplification process is reflected in the linear combinations 
in eqs. || and ^ and is the reason for the rotation of mode polarization. After amplification, any 
polarization the original radiation might have had is irrelevant. For example, even for thermal 



emission in a magnetic field polarized according to eq. 12, each mode contains random linear 
polarization (arbitrary q and u) for propagation in any direction other than 8 = 0. Because the 
projection of every II on the q—u plane points at a random direction, each mode rotates differently 
during the amplification process, individual polarizations are randomized, and the initial overall 
polarization disappears once the radiation field is dominated by the induced photons. The only 
polarization that can survive the amplification process is that of stationary modes that do not 
rotate, locking the individual polarization vectors. All four Stokes parameters then grow at the 
same rate and the fractional polarization remains constant. The stationary modes of amplified 



radiation are the eigenvectors of the matrix K and were identified in E96 . Denote by Kq(x) the 
unsaturated absorption coefficient of the maser transition in the absence of magnetic fields. In the 
presence of a magnetic field, the unsaturated absorption coefficient of the Am transition becomes 

K& m {x) = k {x + x B Am) ~ kq(x) + x B Am dK °} X \ (14) 

ax 

similar to equation [l0| for the corresponding intensities. Again, the last relation is appropriate for 
xb "C 1- Maser polarization properties are controlled by the two ratios 

R = «o +«o _ 1 

2kq Kq(x) dx 

In each case, the last expression is the result to leading order in xb, with the corrections smaller 
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by order x 2 B because of symmetry. The solution for the maser circular polarization when xb <C 1 is 



8R C ^ 8Av B k'q{v) 
(2R~i + 1) cos 9 ~ 3 cos 6»k H' 



(16) 



where 9 is restricted by y/2/3 > cos# > 4\^2xbx (see E96 , eqs. 4.13-4.15). The last expression is 
obtained by inserting R\ and R c from the previous equation, corrections are smaller by order x\. 



4.1. Unsaturated masers 

When the maser is unsaturated (</(f ) <§C </s), to leading order in both J/J s and xg, K m — kq 
and = k c = for any mode, as is evident from equations [5] and 14. In this approximation the 



polarization rotation is neglected and the intensity follows the familiar solution of an unsaturated 
scalar maser I(v) = Sq exp t(z^), where t(v) = kq{i>)1 is the amplification optical depth, irrespective 
of polarization. Therefore I' /I = tk' /kq, and for any polarization v 

= w«M. (17) 

t(v) k q {u) 



Inserting the stationary circular polarization from equation 16, the 7£-profile of the maser solution 
in the unsaturated regime is 

K(v) = f^f (18) 
6 cos 

where 

t[u) , l{y) 

V ; ro + ln-^- 
^(^o) 

The profile ^{v) expresses the spectral shape of (the unmeasurable) t{v) in terms of the measured 
intensity and the free parameter To = t(vq), the maser optical depth at line center. This profile 
increases toward the line wings from a minimum at line center, fundamentally different from the 
flat profile of thermal radiation. When such TZ spectral shape is detected the maser is unsaturated, 
therefore the spectral profile of ln/(z^) can be used to deduce the profile of ^o(^). The kq profile 
determines the intrinsic Doppler linewidth, a quantity which is not known a-priori and which is yet 
to be reliably determined in any maser source. Furthermore, the spectral relation between 1Z{v) 
and \al{u) provides a method for direct determination, the first of its kind, of the maser optical 
depth to from observable quantities. Once tq is found, can be used to determine B/cos9. 

It may be noted that the polarization vector of the maser solution is not fully stationary during 
the early stages of maser growth when J / J s < x\\ indeed, no polarization mode can avoid rotation 
at that stage. An estimate of this early rotation can be obtained by including the terms linear in 
xb while maintaining zeroth order in J/J s - In this approximation order n m remains equal to kq 
so the intensity growth is the same, and k/ still is zero so the linear polarization does not rotate. 
However, k c ^ and the circular polarization does vary according to 

= (l - v 2 ^ R cCO s9. (20) 
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The solution of this equation for any initial polarization v\ is 

Vi + tanh tR c cos 9 

v — 

1 + Vi tanh tR c cos 9 ' 



(21) 



a solution valid as long as exp(To)S'o/J s < x^. Pumping schemes of OH masers typically have 
So/J s ~ 10 -5 and the observed 20% polarization is reproduced with xb = 0.03 for propagation at 
cos 2 9 = 1/3. Therefore this solution holds for To ln(100) = 4.6. Figure 1 displays the 7?.-profiles 
of the polarization solution (v, = vq) at tq = 5 with and without the rotation. The difference 
between the two is practically insignificant. Furthermore, note that the neglected terms of higher 
order in J / J s will push the polarization toward stationary behavior and thus decrease the variation 
of v below the result of equation Therefore, this approximation produces an upper limit on 
variation of v and the actual solution must fall between the two displayed profiles. As tq increases, 
the linear polarization joins the rotation so that the overall rotation of the polarization solution 
decreases further in inverse proportion to J and the stationary approximation becomes even better. 
The deviation from flatness of the maser 1Z profile is a general result, reflecting line narrowing 



during unsaturated growth (e.g., Elitzur 1992 §4.5). It is instructive to examine the unsaturated 
amplification of polarized thermal radiation with an arbitrary intensity profile Ith( l/ ) an d a constant 
7Z profile, i.e., V oc I'. Equal amplification of the Stokes parameters produces 

Uy) = /th(^)e rM , V(y) = Au B cos9lMe T ^. (22) 

Evidently, V can not remain proportional to I' because the derivative of the amplified intensity is 
not equal to the amplified intensity derivative. The proportionality would be retained only if t(u) 
were independent of v instead of being sharply peaked at line center. Since the spectral shape of 
I t h follows the same Doppler profile as Kq, straightforward algebra yields 

Av B cos 6 

1 + Tq exp(-X^) 

The 1Z profile is constant only at tq = 0, i.e., only for the input radiation. The amplification 
process destroys the profile flatness even though it only amplifies the intensity without affecting 
the polarization. When To > 1, 1Z assumes instead a Gaussian absorption shape. This happens 
because the amplification is centered on vq, therefore the / ± components are amplified more strongly 
on their inner shoulders, the ones closer to v§. While this preserves the general symmetry principle 
of eq. 13, as it must, the invariance under frequency shifts of eq. [U] is destroyed because l^{y) are 



no longer symmetric about their centroids at v± = vq ± Avb- 



4.2. Saturated masers 

When the maser is saturated (J{v) > J s ), the radiative transfer equation for the intensity of 
the dominant ray of the polarization solution becomes 

dI (v) o j kq(^) (9A , 

~1T = 37tJs ^ (24) 
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where is the maser beaming angle, which in general may vary with frequency shift from line 
center ( Elitzur 1990| ). Therefore, the profile of I[y) has the spectral shape of the ratio kq(is) /Cl(i>) . 



The frequency variation of depends on the type of maser amplification, which is controlled 
by the geometry (Elitzur, Hollenbach & McKee 1992). In matter-bounded masers, whose prototypes 



are filaments, the beaming angle is independent of frequency and 1(f) oc Kq(u). In amplification- 
bounded masers the beaming angle varies with frequency shift from line center according to dimen- 
sionality. In the prototype planar maser, a saturated disk, £l{y) oc 1/kq(u) and I(v) oc [kq(v)] 2 . And 
in the prototype three-dimensional maser, a saturated sphere, oc 1/[ko(v)] 2 , so I(v) oc [«o( l/ )] 3 - 
In summary, the intensity spectral profile of a saturated maser obeys 

I{u) oc [«o(i/)]f (25) 

where p is the dimensionality of the geometry: 1 for filaments, 2 for planar masers and 3 for 
sphere-like configurations. Therefore I' /I = pk' /kq and the 7£-profile of a saturated maser is 

K{v) = (26) 
3p cos 6 

The unsaturated profile \P is replaced by the constant 1/p during the saturation process, reflecting 
the rebroadening of the maser line. Similar to thermal radiation, 1Z is constant across the saturated 
maser line, only the value of this constant is different. For a given magnetic field and propagation 
at the smallest angle allowed for maser polarization, 9 = cos -1 y/2/3 = 35°, the proportionality 
constant for filamentary masers is 4 times larger than for thermal radiation (twice as large for 
disks, 1.3 times for spheres). In other words, masers require smaller fields to produce the same 
circular polarization as thermal radiation. The disparity between the two cases increases with 9 in 
proportion to 1/ cos 2 9. 

The saturation effect was introduced here in its standard form. Strictly, this form applies only 
to linear masers, three-dimensional effects lead to a more complex form ( Litvak 1973| , Bettwieser 



& Kegel 1974, Neufeld 1992). However, because of the tight beaming of maser radiation, the 
standard form provides an adequate approximation around line center, failing only in the extreme 
wings (x > 2; [Elitzur 1994 ) of saturated masers, and its use is justified in the present analysis. In 



addition, frequency redistribution has been neglected. Incorporation of this important ingredient 
into maser theory is sufficiently difficult that it has not yet been fully accomplished even for scalar 
masers. In keeping with common practice, it has been neglected in this first study. Also, the maser 
polarization problem was solved only for a J = 1 — ► transition. However, when xb <C 1 the 
solution holds for all spins (E96) and thus is applicable here. Note, though, that this result holds 
only for isotropic pumping. Pumping schemes that will introduce m-dependence would require 
separate handling. 



5. RESOLVED ZEEMAN PATTERN; x B > 1 

When xb > 1 the Zeeman components separate and the radiation displays three fully polarized 
lines whose linear and circular fractional polarizations are constant across each component. The 
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circular polarizations of the two a-components, centered on v± = vq ± Avb, are 

v± = ± 2COSS 
1 + cos^ 



whatever the degree of saturation ( Poldreich, Keeley fe Kwan 1973| ; p96| see |Field fc Gray 1994| 



for the approach to this solution of amplified unpolarized background radiation). In fact, this 
result follows from general symmetry properties and holds also for thermal radiation, as can be 
easily verified from expressions listed in Crutcher et al. 1992 . Although the polarization is constant 



since v 




across each individual profile, it maintains the reflection symmetry of eq. 13 
Independent of the nature of the radiation, V(v) is proportional to I{v) across each component, 
not to The ratio profile now obeys 1Z = v^I/I', leading to 

± 2 f 1 thermal 

1Z(v) = —— t — - x \ ^{v) unsaturated maser (28) 

saturated maser 

if kq and the thermal radiation follow a Doppler profile with width Avq. In the profile \& of each 
unsaturated component, the intensity is normalized to its line center magnitude for that component 
(cf. eq. 0). 

The ratio profiles for <C 1 and xb > 1 are fundamentally different from each other. In one 
case 1Z is symmetric in reflections about the center of the single observed line, in the other it is 
anti-symmetric about the center of each observed u-component at v±. This difference provides a 
decisive observational test to determine whether the Zeeman splitting is larger or smaller than the 
linewidth. 



6. OBSERVATIONAL IMPLICATIONS 

Spectral analysis of circular polarization provides a new powerful tool for maser studies. Com- 
parison of the profiles of V{u) and I\v) offers sensitivity to details at the level of ~ xb and has long 
been utilized in studies of thermal radiation. The maser analysis can be readily performed with 
the techniques developed for the thermal case (see e.g. Ireland &: Heiles 1982|) where the observed 
V spectrum is fitted as 

V(u) = al{y) + bl'(u) (29) 

and a and b are free parameters. The al term is introduced to account for differences in instrumental 
sensitivity in left and right circular polarization. The parameter a has no bearing on the source 
properties, it is merely adjusted to make V — al anti-symmetric about line center. The next step 
is to fit V — al to the intensity derivative producing the meaningful parameter b = AubcosO 
(eq.©. 

In the case of maser radiation, various mechanisms can filter one sense of circular polarization 



in the source itself (see e.g. §6.7 of Elitzur 1992 and references therein). Then the observed V profile 



may not be perfectly anti-symmetric even before it reaches the detection instruments. Similar to 
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instrumental effects, the left-right asymmetry introduced by such filters can be handled by fitting 
the observed V profile of an unsaturated maser as 

(30) 



(eq. 18), with a adjusted to make V — al anti-symmetric about line center. Fitting the resulting 
anti-symmetric profile produces the adjustable parameters b = 8Ai/b/3cos# and c = To, the maser 
optical depth at line center. This is the first method to directly determine the maser optical depth 
for self amplification and the intrinsic velocity Doppler width (obtained from the spectral shape of 
In/). 

When the maser is saturated 7Z is flat, so the V profile is fitted as in the thermal case only 
b = 8A^s/3pcos0 (eq. p6| ). The parameter p is determined by the maser geometry (it is 1 for 
filaments, 2 for planar masers and 3 for three dimensional masers), and this cannot be found 
directly from observations. However, strong masers are unlikely to be spherical, so p = 3 can 
usually be dismissed. In particular, the relevant geometry for shock generated masers, applicable 



to water masers as well as the 1720 MHz masers at the Galactic center (Yusef-Zadeh, Uchida & 



Roberts 1995; Frail et al. 1996; Yusef-Zadeh et al. 1996| ), is likely to be planar or filamentary and 



p is either 2 or 1. Another handle on p can come from the linewidth Av m of the observed maser 



intensity. When kq has a Doppler shape with width Ai>d, equation 25 shows that 



Au m = (31) 

and the linewidth observed in disk masers is only 70% of the Doppler width. A line narrower than 
the expected thermal width in a saturated maser can be taken as indication of planar rather than 
filamentary geometry. 

The distinct shape of the unsaturated 1Z profile, evident in figure 1, provides the first direct 
method to determine the saturation stage of a maser. Detection of the predicted logarithmic 
variation across the line would provide a unique, unambiguous signature of unsaturated maser 
operation. However, a flat 1Z cannot immediately be taken as conclusive evidence for saturated 
operation; within the observational errors, such behavior can always be attributed to sufficiently 
large To. If T> denotes the dynamic range of intensity measured across an unsaturated maser line, 
the expected relative variation of 1Z is (l/ro)rnD. The observed profile could be indistinguishable 
from a constant if this variation is smaller than the observational error e, namely, if 

T >iln£>. (32) 

A flat 1Z only implies a lower limit on tq. However, if the maser is indeed unsaturated, then there 
is also an upper limit on tq. In this case, the measured brightness temperature at line center is 
T{, = T x expTQ, where T x is the maser excitation temperature. This quantity in turn is related to 



the fractional inversion r\ = (n2 — ni)/(ri2+ni) viaTj, = huo/2krj (e.g. Elitzur 1992 §4.2). Therefore 



tq = \a{2r]kTi ) / huo) and an absolute upper limit on tq is obtained by inserting in this expression 
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rj = 1, the maximum theoretical value, yielding 

2kT b , . 

r <ln— A 33 

The optical depth of an unsaturated maser is always bounded by the observed brightness tempera- 
ture. The maser can display a flat 1Z profile only if this upper limit is compatible with the previous 
lower limit. 

To summarize the observation analysis procedure: The free parameter a is varied until the 
combination V — al of the observed V and I profiles is anti-symmetric about line center, then the 
ratio I'/ (V — al) is formed. If this ratio displays linear variation with lnl across the line then the 
maser is unsaturated and the slope and intercept of the linear fit determine B/cosO and tq. If, 
instead, the ratio I'/(V — al) is constant across the line then the situation depends on the last two 
bounds. If these bounds are incompatible, the maser is saturated. That is, if the fitting error e and 
the observed brightness temperature Tt> and dynamic range T> obey 

ln2?>eln^ (34) 

yet I'/{V — al) is constant across the line then the maser must be saturated. When the reverse 
relation exists between T>, e and T&, the saturation stage is unsettled; the data lack the sensitivity 
to distinguish between the 7Z profiles of saturated and unsaturated maser operation. 

Observational data that obey the last inequality are of sufficient quality to determine conclu- 
sively whether the maser is saturated or not. The recent VLA polarization data of OH 1720 MHz 
masers near the Galactic center by Yuscf-Zadeh et al. 1996| are the first to meet this criterion. The 



brightness temperatures of all resolved Galactic OH masers (hv^jk = 0.08 K) always obey Tf, ;> 10 10 
K. If any such maser is unsaturated, the upper bound of eq. ^ implies that tq < 26 (a realistic 
estimate for the inversion, 7] <^ 10%, would reduce this bound to to < 24). The 1Z profiles of the 
Galactic-center masers were fitted as constants by Yusef-Zadeh et al. 1996| and the errors listed for 



the magnetic fields are less than 10% in all sources, as low as 3% in the strongest (source A). These 
are essentially error estimates of the quality of the fits to flat ratio profiles, thus e ~ 0.03-0.1. The 
resulting fits match the observations all the way to the noise level of ~ 10 mJy for an intensity 
dynamic range T> ^ 1000 in the strongest source and no less than ~ 20 in the weakest (source 



F). Therefore, from equation 32 the lower limits on To range from 30 to 230 and are incompatible 
with the upper limit for unsaturated OH masers. This conclusion can be strengthened by a proper 
analysis that would test the hypothesis of an unsaturated TZ spectral shape with the observations. 
Such analysis should produce constraints that are even tighter and more secure. Still, it seems safe 
to conclude already on the basis of the current analysis that at least the strongest masers, and quite 
possibly all of them, are saturated. Previously, the strongest evidence for saturation was provided 



by OH 1612 MHz masers in late-type stars ( Elitzur 1992 , §8.6). This maser emission follows the 



temporal variation of the IR radiation, the presumed pump, and the similar amplitudes indicate a 
linear relation between pump and maser intensity, as expected in saturated operation. While this 
evidence seems persuasive, it is indirect as it requires an assumption, however plausible, regarding 
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the pump mechanism. In contrast, the new evidence is direct, obtained from the relation between 
I and V, two maser intensities measured independently. 

Detailed profile analysis of circular polarization also enables direct, unambiguous determina- 
tion of the relation between the magnitudes of the Zeeman splitting and the line Doppler width. 
Substantial circular polarization (50% and higher) has long been observed in OH maser emission 
from late-type stars (e.g., [Reid et al. 197S , Claussen fc Fix 1982 , Cohen et al. 1987 ). In the absence 



of theory for xb < 1 at the time, this was taken as a signature of xb <; 1, implying magnetic fields 
of at least ~ 1-10 milligauss. However, this polarization often appears as sharp reversals between 
adjacent narrow spectral components of the parameter V, as expected for xb < 1- Analysis in 



terms of the theory for xb < 1 yields instead fields of only ~ 0.1 milligauss (E96), so it is impor- 
tant to determine conclusively what is the relevant domain of xb ■ The new profile analysis provides 
a simple, reliable method to resolve this issue. If the V profile can be fitted across each component 
as V = al without residuals then xb > 1. If, on the other hand, the residuals are significant and 
V — al is anti-symmetric around line center then xb < 1- 

Finally, circular polarization analysis, of either thermal or maser radiation, can never fully 
determine the magnetic field because the propagation angle 9 is not known. However, in the case 
of maser radiation 9 can be determined from linear polarization measurements; the polarization 
solution shows that both circular and linear polarizations are generated when xb <C 1, and to 
leading order in xb the linear polarization is a unique function of 9 ( |E96| , eq. 4.14). Indeed, strong 
linear polarization has been recently reported for the 1720 MHz masers by |Clausscn et al. 1997 . 



Although Faraday rotation can reduce the linear polarization, even eliminate it on occasion, study 
of linear polarization is an important task for observations. 

Discussions with Tom Troland and Farhad Yusef-Zadeh, and the partial support of NASA 
grants NAG 5-3010 and NAG 5-7031 are gratefully acknowledged. 
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Fig. 1. — The 1Z profile of the polarization solution for an unsaturated maser with the parameters 
listed. The full line is the profile including the maximum possible polarization rotation (eq. 21), 
the dashed line is the profile neglecting this rotation (eq. 19). 



